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Abstract. We consider classes of simply connected planar domains which 
are isophasal, ie, have the same scattering phase s(A) for all A > 0. This is a 
scattering-theoretic analogue of isospectral domains. Using the heat invariants 
and the determinant of the Laplacian, Osgood, Phillips and Sarnak showed 
that each isospectral class is sequentially compact in a natural C°° topology. 
This followed earlier work of Melrose who showed that the set of curvature 
functions k(s) is compact in C°°. 

In this paper, we show sequential compactness of each isophasal class of 
domains. To do this we define the determinant of the exterior Laplacian and 
use it together with the heat invariants (the heat invariants and the determi- 
nant being isophasal invariants). We show that the determinant of the interior 
and exterior Laplacians satisfy a Burghelea-Friedlander-Kappeler type surgery 
formula. This allows a reduction to a problem on bounded domains for which 
the methods of Osgood, Phillips and Sarnak can be adapted. 



1. Introduction 

1.1. The isospectral problem. In this paper, we consider a scattering-theoretic 
version of the famous question 'Can one hear the shape of a drum?' posed by 
M. Kac pT0| . In mathematical terms the question is whether a planar domain O 
is determined up to isometry by its Laplace spectrum (with Dirichlet boundary 
conditions, say), where the spectrum < < A| . . . is counted with multiplicity. 
The answer to this question is known to be negative || (though there are some 
positive results for restricted classes of domains P^l). In view of this, it is reasonable 
to ask how 'small' is the set of domains isospectral to a given domain. One way to 
make this precise is to ask whether the isospectral class is compact in some topology 
on domains. Melrose |Q showed that this is the case, where the topology is taken as 
the C°° topology on the curvature function k(s) : s G [0, L] of the boundary of the 
domain [L is fixed over the isospectral class, as discussed below). A disadvantage of 
this topology is that it does not exclude the possibility of a sequence of isospectral 



domains pinching off (see figure 1.1) 



This result was proved using the 'heat invariants'. The heat invariants of a 
domain O are coefficients in an expansion of the heat trace ee>(i) as t — ► 0. Since 



e (i) = tre-* Ao = £\ 
4=1 
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FIGURE 1.1. Pinching off without blowup of curvature 



the heat trace is a spectral invariant. It has a well known asymptotic expansion 

oo 

(1.1) eo(t)~ a j t ~ j/2 > 

J=-2 

as t tends to zero, with the heat invariants Oj 'local', that is, integrals over the 
domain O or the boundary H = dO of locally defined geometric quantities. The 
first few are 

area(O) 1 f 

a - 2 -^r~ = ' 1 



(1.2) 



_i = c_i length(<90) = c_i / Ids 

J do 

ao = cox(O) = c / fc(s) ds 



a\ = c\ I k 2 (s) ds 

JdO 

with Ci 7^ 0. Melrose showed that 



(1.3) a 2 i-i = Cl (k (l \s)) 2 + Pl (k(s),...kV- 1 \s))ds ci^O 

JdO 

where pi is a polynomial; this is the main step in his result. 

Notice that the first two heat invariants give A = area(O) and L — length(<90) 
as spectral invariants. Consideration of the isoperimetric quotient shows that the 
disc of radius r, D r , is determined by its spectrum, an observation that perhaps 
led to Kac's question. 

Osgood, Phillips and Sarnak (abbreviated OPS from here on) considered the 
question of C°° compactness from a different point of view Jl6| , Jl7| . They regarded 
a planar domain as the image of the unit disc D under a conformal map F . The 
metric on the domain is then isometric to e 2 ^go where go is the flat metric on the 
disc, and <p — log \F'\ is a harmonic function. Thus <f> is determined by its boundary 
values 4> \ dD. Given a harmonic function 0, one can find the corresponding domain 
F(D), which is a flat planar domain (possibly self-overlapping). OPS showed that 
isospectral classes are compact in the C°° topology of <fi restricted to dD. This 
result excludes degenerations of the form illustrated in figure |l.l| since derivatives 
of F must blow up under such a degeneration. (Melrose also has an argument ruling 
out such degenerations using the first positive singularity of the wave trace 

Osgood, Phillips and Sarnak used the heat invariants, plus one other invariant, 
the determinant of the Laplacian (described below), to deduce their result. They 
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exploit a formula, due to Polyakov (T3| and Alvarez |2| , expressing the determinant 
in terms of the function <f>: 

(1.4) logdetAo = -L / <t>d n <t>-^- f <j> + logdet A D . 

12tt J s i 6tt J S i 

The formula is remarkable since the first term on the right is nonnegative and 
almost the square of the Sobolev \ norm of <f>. 

1.2. An analogous problem for the exterior domain. In this paper we are 
interested in the exterior Laplacian. Let = M 2 \C, the exterior of an obstacle, and 
let An be the Laplacian on L 2 (f2) with domain H 2 (fl) n H^fl) (ie, the Dirichlet 
Laplacian). The operator Aq is self-adjoint with continuous spectrum on [0, oo) 
( pijl , chapter 8). We wish to formulate a problem about exterior domains that is 
analogous to the isospectral problem. 

To do this, we observe that the isospectral condition may be expressed in terms 
of the counting function 

Ao(A) = number of eigenvalues of Aq < A 2 . 

To say that two domains have the same spectrum, counted with multiplicity, is 
equivalent to saying that they have the same counting function N(\), and the 
problem considered by OPS is to show that the class of domains with a fixed N(X) 
are compact in some natural topology. 

For exterior domains, = R 2 \ 0, it is known that 1/2-7T times the scatter- 
ing phase s(A) is analogous to the counting function. The usual definition of the 
scattering phase is 

s(A) = —i log det S'o(A), 
where S'o(A) is the scattering matrix (see [^5|). For our purposes, it is more illu- 
minating to note that the difference between the spectral projection Eq(\) on the 
interval (— oo, A) for An, and the corresponding spectral projection Eo(X) for Ar2, 
is trace class in a distributional sense, with 

(L5) oo 

tr f <f>'{<T)(En(<T)-E (a))d(T = tr(<f>{An)®0-<f>(Aw))= f <j>' {<r)^^- da 
Jo Jo Zn 

for any </> S C£°(R) (see ||; the normalization of their scattering phase 6(X) is 

minus one-half of our s(A) — cf remark 1 of their introduction). Thus — s(A)/27r is 

a regularized trace of the spectral measure. Since 

No(Vo) =trE a (a), 

the analogy between the counting function and the scattering phase is clear. Let 
us say that two obstacles are isophasal if they have the same scattering phase. 

A strong indication that it might be possible to use the scattering phase to prove 
compactness results about isophasal classes of domains comes from noting that the 



formula (1.5) holds also for (j)(a) = e . Thus the regularized trace of the heat 



kernel is given in terms of the scattering phase by 
(1.6) r-tre-* An = tr ( e - tAn - e -*<V) = -- / s(X)e~ x ^ f \d\. 



However, direct construction of a parametrix for the heat kernel of Aq near t = 
shows that the regularized trace has an asymptotic expansion of the form ( |1.1| ) with 
the same coefficients (up to changes of sign). Thus we immediately get Melrose's 
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result for C°° compactness of the curvature function of the boundary. The question 
is then whether this can be improved, OPS-style, to a result of C°° compactness 
of the domain. It is very natural to look for an analogue of the determinant of the 
exterior operator in order to do this. 

1.3. Determinants and surgery formulae. We begin by recalling the definition 
of the determinant. Let A be a strictly positive elliptic mth order differential 
operator on a bounded domain of dimension n (compact manifold, possibly with 
boundary). Then A has positive, discrete spectrum < fj,\ < fj,2 • ■ • — > oo. The 
determinant of A is defined in terms of the zeta function, £(s). The zeta function 
is defined by 

M C<»> = £^ = r^f 

in the region of absolute convergence, Sfts > n/m. Since the heat trace has an 
expansion pQj 

oo 

(1.8) tre-* A ~ t-^ m bj, t->0, 

j=—n 

it follows that J t s_1 tr e~ tA dt continues meromorphically to the complex plane 
with at most simple poles at s = —j/m,j > —n. The factor r(s) -1 vanishes at 
s = ensuring that the zeta function is regular at s = 0. The determinant of A is 
then defined by 

logdet^ = -C'(0). 

If A is not strictly positive, that is, has a zero eigenvalue, then the determinant is 
defined to be zero. However, it is usually of interest to look instead at the modified 
determinant, det 'A. This is defined by defining the zeta function using only the 
nonzero eigenvalues of A, and then taking log det 'A — — £'(0). An equivalent 
definition is that the modified determinant of A is the determinant of A + Ho, 
where IIo is orthogonal projection onto the null space of A. 

If A is pseudodifferential, then the definition above does not make sense in gen- 
eral, since the heat trace of A may have log terms (terms of the form t^ m \ogt) 
as t — > 0, and then the zeta function may have a pole at s = 0. However, in the 



case of interest in this paper — the Neumann jump operator (see Definition 2.5) — 
one can rule this out and then the log determinant is defined just as for differential 
operators. 

The terminology 'determinant' is justified by the fact that if A were an operator 
on a finite dimensional space, and therefore had a finite number of eigenvalues, then 
we would have 

-C'(0) = - X)( _1 °8/*j)A*7"l«=o = X)l°gA*j = log det A 



The log determinant is a non-local quantity; that is, it cannot be written as 
the integral over O or dO of locally-defined geometric quantities JiiJ. However, 
it behaves in many situations as a 'quasi-local' quantity, in the following sense: 
when a localized perturbation is made in the operator, one can often find a formula 
for the change of the log determinant which involves only the perturbation. An 
example, which is highly relevant to this paper, is the Mayer- Vietoris type surgery 
formula for the log determinant proved by Burghelea, Friedlander and Kappeler [0 
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(henceforth BFK). Suppose that A is an elliptic partial differential operator on a 
compact manifold M, and H is a hypersurface, with M the manifold with boundary 
obtained by cutting M at H . Let B be an elliptic boundary condition for A on the 
boundary of M. BFK found a formula for logdet A — logdet(A, B) in terms of the 
log det of a pseudodifferential operator R on H and other data defined on H . In 
the particular case of the Laplacian on a two-dimensional manifold M, with H a 
curve dividing M into two components M\ and M2, and B the Dirichlct boundary 
condition, they showed that 

(1.9) log det A M - log det(A Ml , B) - log det(A M2 , B) = log det R + log a - log /, 



where R is the Neumann jump operator (see Definition 2.5), a is the area of M and 
I is the length of H . 

In this paper, we look at the exterior and interior Dirichlet Laplacians for an 
obstacle O from this point of view. Thus, we consider the boundary H of O to be 
a cutting of the manifold R 2 and look for a BFK-type surgery formula. Of course, 
one problem is that R 2 is unbounded, so An has continuous spectrum and its log 
determinant is not defined. This is the topic of the next section. The main theorem 



is Theorem 2J3, wh ich gives a surgery formula for this regularized log determinant 



very similar to fll.9| ). The proof of this theorem is the subject of the third section. 
1.4. Compactness of isophasal sets. In the fourth section we show that each 



class of isophasal sets is compact in a natural C°° topology (Theorem 4.1). First 
we must specify the topology on domains. Following Osgood, Phillips and Sarnak, 
we define a sequential topology, ie, we specify the convergent sequences rather than 
the open sets. This is appropriate since our goal is to prove sequential compactness. 

It is inconvenient to deal with unbounded domains, so we pass to the inversion 
f2 7 of f2. We will say that a sequence of exterior domains tti converges in the C°° 
topology if there are Euclidean motions Ek of the plane such that 

(i) the closure of Ek^lk does not contain the origin; 

(ii) the sequence Ok of inversions of Sfcfife about the origin converges in the sense 
that there are conformal maps Fk from the disc to Ok-, with \FL\ never zero, which 
converge in the H s topology for all s. 

The notion of convergence of Ok in condition (ii) is stronger than convergence 
in the OPS topology, which would say that there are Euclidean motions Ek and 
conformal maps Fk from the disc to EkOk that converge in H s for all s. It is 
important that the group of Euclidean motions is allowed to act on the domains 
VLk and not on the Ok- 



To prove Theorem 1.1, we use the conformal equivariance of the Laplacian in 
two dimensions to compactify the problem. That is, we consider a metric g on R 2 
which is a conformal multiple of the flat metric, so that infinity is compactified to 
a point. Comparing our surgery formula to that of BFK on the compactified space, 
we show that the Laplacian on a certain bounded domain of S* 2 (depending on the 
obstacle) has fixed determinant, as the obstacle ranges over an isophasal set. This 
allows us to adapt the argument of OPS to obtain the result. 



2. DETERMINANT OF THE EXTERIOR OPERATOR 



In this section we shall define the modified determinant of the exterior Laplacian, 
and state the main theorem. 
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The determinant is usually defined in terms of the zeta function, which in turn 
is defined in terms of the trace of the heat kernel. In the case of the exterior Lapla- 
cian, the heat kernel is certainly not trace class, since it has continuous spectrum. 
However, as discussed in section [L^, the difference between the exterior heat op- 
erator and the free heat operator is trace class for every t (see || and ^); we will 
denote this trace by r-tre~ tAn , and call it the regularized heat trace. Thus, the 
obvious candidate for the zeta function is 

(2.1) feW — ^jff^.— f 

Unfortunately, this integral does not converge for any value of s. To deal with 
this we break up the regularized heat trace into two pieces. Let ^ be a smooth 
f unc tion that it identically one near A = and identically zero for A > 1. Then, by 
(U), r-tre~* Af! = ex(t) +e 2 (t), where 

(2.2) e 1 (t)=— / X (A) S (A)e^ 2 *AdA 

7T 







and 



t f°° 

(2.3) ea(t) = / (1 - x(A))s(A) e - A2 *AdA. 

t Jo 

We write £n(s) — Cn , l ( s ) + Cfi2( s ) f° r the corresponding decomposition of the zeta 
function. Then e2{t) is exponentially decreasing at infinity. On the other hand, 
since r-tr has the usual asymptotic expansion at t = 0, and e\{t) is smooth at 
t = 0, we see that e2(t) has an expansion of the form (fL~8|). Thus, 

(2.4) ^ = fw/o tSC2{t) - 

continues meromorphically to the entire plane with no pole at s = by the usual 
argument. The other part, Cn,i( s ) rnay be directly expressed in terms of the scat- 
tering phase by the formula 

(2-5) ( nA ( s )=— A- 2 MA) X (A)-. 

k Jo A 

To understand this integral we need the asymptotics of s(A) as A — 0. Following 
|6j, let us define, for this and other purposes, the function ilg. 

Definition 2.1. The function ilg A is defined to be 

log(l/A) 

it goes to zero as A — > 0, but slower than any positive power of A. 

Lemma 2.2. For any obstacle O, the scattering phase satisfies 
(2.6) s(A) =vrilgA + 0((ilgA) 2 ), A - 0. 

Here, the 0((ilgA) 2 ) term is uniform over each isophasal class. 
Proof. In the appendix we compute s(A) for a disc of radius 1; the result is 
s Dl (X) =^il g A + 0((ilgA) 2 ), A-0. 



DETERMINANTS OF LAPLACIANS IN EXTERIOR DOMAINS 



7 



The scattering phase for a disc of radius r is then sjj r (X) = srj 1 (r 2 X). Thus, for a 
disc of radius r, we also have 

| SDr (A)-7rilgA| <C(r)(ilgA) 2 . 

Using the first two heat invariants, the area and perimeter are constant on an 



isophasal class of domains; hence, by Lemma 5.1 the inradius and circumradius are 



uniformly bounded below and above. Thus, we can sandwich any domain in an 
isophasal class between fixed discs D r and Dr. By s(A) is monotonic in the 
domain, so we obtain (2.6). □ 



Substituting this expansion into (2.5), we see that the zeta function is meromor- 
phic in the half plane 5is < 0, but not in any neighbourhood of s = 0. To see this, 
consider the function 

roc j\ 

9(s) = J A- 2 *ilgA x (A)-. 

By differentiating once in s, it is not hard to show that g(s) is equal to — log(— s) 
plus a smooth function as s f 0. Thus, the zeta function has an expansion of the 
form 

(2.7) Cfi(s) = ao- s\og(-s) + a 2 s + 0{s 2 logs), s t 0. 

This allows us to make 

Definition 2.3. The logarithm of the determinant of the exterior Laplacian is 
defined to be 

log det 'Ao = -a 2 , 



where a 2 is the coefficient of s in the expansion (2.7) 



For future use, we observe here that if we consider the operator Ao + \i instead 
of Ao, with fi > 0, then the zeta function is given by 

r°° rll 
(2-8) CM*) ' 



t s r -tre- tAn e^*— ; 

t ' 



the integral is now defined for 5Rs > 1 and continues meromorphically to the complex 
plane with no pole at s = 0, since the exponential factor e - ' 1 ' makes the integral 
convergent at infinity for any s. However, it is useful to write the zeta function in 
the same way as for the case /x = 0: 



(2-9) CnA*) 



(A 2 + /i)- s - 1 x (A)s(A)AdA + 



The determinant is then defined in the usual way. 

Definition 2.4. For fj, > 0, the logarithm of the determinant of Ao + M is defined 

by 

logdet(Ao+^) = -CM0). 

Remark. We write det ' instead of det in definition |2.3| because it is more similar 
to the modified determinant described in section 1.3 than the determinant — this 



becomes clear in the calculation of section 3.3 



It is not yet clear that the quantity in Definition 2.3 merits the term 'determi- 
nant'. We believe that the following theorem justifies the definition — compare with 



equation (1.9). First we give a formal definition of the Neumann jump operator R. 
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Before stating it, we observe that for any jj, > 0, and given any continuous function 
/ on H , there is a unique bounded extension u of / to Q satisfying (An + fi)u = 0. 

Definition 2.5. The Neumann jump operator R for the obstacle O is the operator 

/ i-> d v u x - d v u 2l 

where / 6 C Q (H), ui, respectively 1*2 are the bounded extensions of / to O, 
respectively f2 satisfying Aui = 0, and v is the outward normal. (This choice of 
normal means that R is a nonnegative operator.) The operator R(fJ-) is defined 
similarly, replacing A with A + fj,. The operators R and R(fi) are pseudodifferential 
operators of order 1, and for /1 > 0, R{fi) is strictly positive. 

Theorem 2.6. The following formula holds: 

(2.10) logdet'An + logdet A a +logdet'i? = 7 + log- 

where 7 is Euler's constant. 

The proof of this theorem is the subject of the next section. 



3. Proof of the surgery formula 



In this section we prove Theorem 2.6. We follow closely the scheme of BFK's 
proof. Thus, the proof consists of three steps. The first step is to establish the 
variational formula 



(3.1) 



— (logdet(A 
djj, 



fi) + logdet(A + fi) + logdet.R0u)) = 0, 



for fj, > 0, where R(fi) is the Neumann jump operator. This calculation was first 
done by Forman Q and the proof given here is almost identical, but it is written 
out in full for the reader's convenience. Thus, integrating (3.1), we find that 



(3.2) 



logdet(A n + M)+logdet(A + fx) + log det H(/z) = C. 



In the second step, we show that C = 0. To do this, we send [i to infinity. 
Then each of the log determinants has an asymptotic expansion in fj,, with local 
coefficients. Clear ly the coefficients of each term must agree on the left and right 
hand side of (3^2) so if we know the coefficient of the constant term for each log 
determinant, then we deduce the value of C . It turns out that the constant term 
in the expansion for each log determinant is zero, so C = 0. 

The third step is to consider the limit \i —* 0. Here we prove the following 
asymptotic expansions: 



(3.3) logdet(A + n) = loglog^ 1/2 + logdet 'A n 

(3.4) log det(A c + fi) = logdet A a +o(l), fi—>0; 

logdet R((i) = -loglog^~ 1/2 



(3.5) 



7 -log2 + o(l), awO; 



L 



log det 'i? - log — + o(l), ^ ^ 0. 
2n 



It is easy to see that Theorem 2.6 follows from this. In the rest of this section we 
give the details of the proof. 
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3.1. Variational formula. First we give a couple of lemmas which will help to 
establish the result. 

Lemma 3.1. For fi > ; the operator (A® + /i) -1 — (A R 2 + /i) -1 is trace class, the 
derivative of 

logdet(A© 

with respect to fi exists, and 

(3.6) -j- logdet(A e +fi) = tr ((A© + /i)" 1 - (A R2 + m)" 1 ). 

Proof. See appendix. □ 

Next we need to introduce some notation. We define the Dirichlet and trans- 
mission Poisson operators, Pdw{n) and Pt r (fi), mapping from H 3 / 2 (H) to H 2 {0)® 
H 2 {n), respectively H^ 2 (H) to H 2 (0) ® H 2 (n) by 

-Pdir(A*)(/) = u i where (A® + jj)u — and u \ H = f 
^ ^ P tr (fi)(f) — u, where (A® + fj,)u = 0, u is continuous at H, [d v ii] = /. 

Here [] denotes the jump in the argument at H (the sign is specified in Definition |2.5[ 
These are the Poisson operators (more precisely, 'half of the Poisson operators) for 
the Dirichlet and transmission boundary conditions considered by Forman. Notice 
that both map to the space 

(3.8) {u E H 2 (0) 8 H 2 {VL) | u is continuous at H}. 



Also note that both (A R 2+/i) * and (A®+/i) 1 map L 2 to (3.8). The corresponding 
trace operators, defined on (3.8), are 

T dil (u) = u\ H 

(3 ' 9) T tT (u) = [d v u]. 

Thus, R(li) = T tI P dil (n) and R(n)- 1 = T dir P tr ( M ). 
Then the following relations hold. 

Lemma 3.2. For [i > 0, 

(3.10) = -(A e + fi^PM, 
a/j, 

(3.11) Pdh(M)T dir P tr (M) = PM, 
and 

(3.12) P tr ( M )T tr (A© + nV 1 = (A© + m)- 1 - (A R 2 + m)- 1 . 



Proof. These are all routine. To prove (3.10), let Pdir (/•*)/ = u{n) and differentiate 
the equations 

(A© + = 0, u{p)\H = f 

and let (d/dfi)u = v to get 

(A© = -i»(m), v(/i)f# = 0. 

It follows that 

(^Piir)/ = = -(A© + = "(A© + ri^PMf, 

which establishes ( |3.1t]| ). 
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To prove the next formula, let u = P tI (j-i)f. Then u is continuous at H; let 
.9 = Idir«. Then Pdi r (£t)<7 is the unique solution v of (A® + fx)v = which is 



continuous at H and takes the value g there. Thus v — u. This demonstates (3.11). 

To prove the final equation, consider the right hand side applied to w. This gives 
us a function u such that (A® +fi)u = 0, with [d v u] equal to T tr (A® + [i)~^w (since 



T tr (A R2 + ii)- x w = 0). Therefore, u = P tr (T tr (A ffi + proving ( p2| ). □ 



Finally we recall from [|] that (d/d/j,)R(fx) is a pseudodifferential operator of 
order —1 for /i > 0, so R(n)~ 1 (d/dfi)R(fi) is an operator of order —2 and hence of 
trace class. Thus, we can calculate 

- j- ( log det flfc)) = - tr (i?( M )-i |-i?( M )) 

= - tr (T dir P tr (/i)T tr ^-P dir (/x)) 

= tr (T dir P tr (M)T tr (A e + /i)- 1 P dir (/i)) by Q 
= tr (P dir ( M )T dir P tr ( M )T tr (A e + M)^ 1 ) 



tr P tr ( M )T tr (A e + M) _i by (3.11) 



tr((A ffi + M )- 1 -(A R2+Ai )- 1 ) In- i-.U-l) 

dfx ' 



A(iogdet(A e +/i)) by (|D 



This completes Step 1. 



3.2. Asymptotics as /i tends to infinity. In this step we calculate the constant 
term in the asymptotic expansion of the log determinants of Aq + fj,, An + /j, and 
R(n) as /z — > oo. In fact, the interior Laplacian and the Neumann jump operator 
have been treated in BFK, where it is shown that in both cases the constant term 
is zero, so we only need to deal with the exterior Laplacian. We use the formula for 
the zeta function in terms of the regularized heat trace to deduce the result. This 
method does not generalize very far, since it requires that the dependence on /i is 
of the form A + /i, but it has the advantage of being very explicit. 

Proposition 3.3. The logarithm of the determinant of An + /! has an expansion 

oo 
J=-2 

as fi — > oo, with po = 0. 

Proof. Recalling (|2.8| ), the zeta function for An + /j is 
(3-13) CMs) = _L/ fr-tre-^e^f. 

r(s) Jo t 

This equation shows why there is an expansion as (i — > oo with local coefficients: the 
factor of e~ M * means that the integral from a to infinity is exponentially decreasing 
in /X, for any a > 0, so only the expansion of the regularized heat trace at t — will 
contribute to polynomial-order asymptotics in /i, and this expansion is local. 



DETERMINANTS OF LAPLACIANS IN EXTERIOR DOMAINS 



11 



For precisely, for any integer k we consider the expansion to 2k + 2 terms of the 
regularized heat trace (see (|L1|)), 



(3.14) 



2k-l 

£< 

3=-2 



at t = 0. Let e&(i) be the difference between the regularized heat trace of e~* Af! 
and this finite expansion. Then, efc(i) is 0(t k ) as i — > and, since the heat trace 
is bounded as t — > oo, efc(i) is also 0(t k ) at infinity. It is easy to see that if eu is 
substituted for the regularized heat trace in (3.13) then both the result, and the 
derivative in s of the result, is 0(ii~ k ) as [i — > oo. Thus, to compute the expansion 
as /i — > oo to this order we need only substitute ( |3. 14 ) in to (3.13), namely 



1 



r( s ) 



aj t j/2 e 



fit 



dt 



Changing variable to t = i/x, this gives us 

r( s + j/2) 



T(s) 



-J/2 



Differentiating at s = gives us an expansion of the form above, with pq = since 
the two T- factors cancel when j = to give a constant. 

□ 



3.3. Expansion as fi tends to zero. Here we consider the asymptotic expansion 
of the log determinants ( |3.3[ ) — (3.5). The second of these, logdet(Ae>+/i), is simply 
continuous as fi — > 0, since Ap + fi has discrete spectrum uniformly bounded away 
from zero as \i — > 0. Thus (3.4) is obvious. 

To understand the behaviour of logdet(An + //), recall from (2.9) that the zeta 
function is given by the analytic continuation of 



(3.15) Cfi,^) 



(a 2 + [iy s - 1 x(\)s(\)\d\ + 



1 



r( s ) 



t s e 2 (t)e-^ 



dt 
T 



from 5fts small, respectively 3?s large. The contribution to the log determinant from 
the second piece is continuous in /i as /i — > 0, so we get precisely —(' n 2 (0) (see 
(U)) in the limit. 

In the first piece, the integrand is compactly supported, and convergent uniformly 
near s = for fixed /.i > 0, since we have the estimate (^]|). So the contribution to 
the log determinant is equal to 



(3.16) 



1 



A 2 



A2 - S (A) X (A)^. 

H A 



Let us write s(A) = 7rilgA + s(A), where, by L emm a 2.2 , s(A) is Q((ilgA) ) as 
A — > 0. (The function ilg is defined in Definition 2.1.) Replacing s by s in ( 3.16 ) 
makes the integral convergent uniformly down to /i = 0, and we get a contribution 
of 



(3.17) 
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}' 



-2 



— ilgA]g = O(ilg^). 



It remains to consider what happens when s is replaced by ir ilg A. Thus, we are 
interested in the expansion of the integral (where for convenience we replace /i by 

f°° A 2 d\ 
(3.18) I ^il S Ax(A) T . 

To calculate this, we break up the integral into pieces. First consider the integral 
from to v. We can estimate the absolute value by 

4ilgAx(A)dA = v-" 1 [ AilgAdA. 
v Jo 

Since (A 2 ilg A/2)' > A ilg A on the interval [0, v] : for small v, this is estimated by 

r A 2 
2 

Thus this term is o{\) as v — > 0, and can be ignored. 

Next consider the integral from v to infinity of ( 3.18| ). We claim that, up to 
an 0(i\gv) error, we can replace the factor A 2 (A 2 + v 1 )^ 1 by 1. To see this, we 
estimate the difference 

,,2 l\ 1-2 

ilg A X (A) - < v 2 j A" 3 ilg A dX. 

Observe that, for small S, A~ 3 ilg A < A -3 ilg A(2 — ilg A) on [0, S], and the quantity 
on the right hand side is equal to the derivative of —A -2 ilg A. Therefore, this term 
is estimated by 

v 2 [(-A~ 2 ilgA)] 2 =0(il gi ,). 
Hence, up to o(l) errors we are left with 



J ilg X— + J ilgAx(A) — . 



Let a = ilg A, e = ilg 5 and x(a) = x(A). In these variables we have 



da f 00 , „ <ia 
— + / *(«)— 

ilg v J e 

= -logilgi/+ ( / x(a)— -log-) 

/•OC J 

= loglog M - 1 / 2 +HR-/ x(«) — , 



(3.19) 



where the last integral is a Hadamard regularized integral (see appendix). This 
may be combined with (3.17) to give 

(3.20) 

f°° do, 

logdet(A n + M ) - loglog M - 1 /2 + HRjf _ c ^( 0) + Q + (i). 

We need to compare this to the log determinant of Aq. By definition, this is 
the coefficient of — s in the expansion of the zeta function as s — ► 0. Recall that 



the zeta function is equal to Cn,i(s) + Cn,2(s) as in fl2.5| ) and ( p.4[ ). Note that the 
contribution from f q 2 is just — 2(^)1 matching one of the terms in the expansion 
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for Aq + From £n,i we get a contribution which is the constant term, as s — > 0, 
of 



A^ s S (A)x(A) 



dX 

T' 



Writing s(A) = 7rilgA + s(A) as before, with s replaced by s above, the inte gral i s 
convergent uniformly down to s = and we get a contribution of exactly ( |3.17 ). 
Thus, it remains to find the constant term in the expansion of 



(3.21) 



\- 2s ilgAx(A) 



dX 



Let us write r = — s, so r > 0. Let a — ilgA; then A 2s = e 2r l a and da/a — 
i\gXdX/X. Let e > be arbitary. Then the integral ( 3.2l| ) is the same as 



-2r/a~ 



da 



da 



r/e 



r/e 



-2r/a~ 



da 



For small e, the factor x(a) may be replaced by one in the last two integrals. 
Writing (3 = r/a, we get 



,-2r/a~ 



r/e 



-2/3^ 

/3 



-2/3 



d/3 



It is clear that there is a divergent term — logr as r — > 0. We seek the limit when 
this divergent term is subtracted from (3.21). This limit is equal to 

f + 2 log, - jfV " - 1)| + jT 00 e -v| 

for every e > 0. Taking the limit as e — > 0, the third term disappears and one 
factor of log e combines with each of the integrals to give two Hadamard-regularized 
integrals. Therefore, we have shown that (3.21) has an expansion 



(3.22) 



log r + HR 



- 2 ^+HR- 



. i do? 

X (a — +o 1 , r^0. 

a 



The first regularized integral appearing here is equal to 7 + log 2 where 7 is Euler's 
constant (pL chapter 5). The seco nd re gularized integral is the same one that 
appeared earlier. Combining it with ( 3.17|) , we get the formula 

logdet'An = 7 + log 2 + HR-/ Cn,a(0) + (@- 



Comparing this with (3.20), we obtain (|3.3|). 



Next we show (3.5). We follow the method of BFK. Since exactly one eigenvalue, 
say Ao(/x), approaches zero as /x — > 0, we have 

logdeti?(/x) = logA (/i) + logdet'i? + o(l) 

as /1 — > 0. Thus, we need to find the expansion of Ao(/x); we use BFK's characteri- 
zation that 

Ao(m) -1 = operator norm of R((j,)~ . 
The operator i?(i/ 2 ) -1 is given by 

w 1 ^ T dil P tI (v 2 )uj = T dir (A M 



14 



ANDREW HASSELL AND STEVE ZELDITCH 



Here J is the map u> <— > uj5h, where Sh is the delta function supported on H . Let 
us analyze the behaviour as v — > 0. By direct computation, it is not hard to show 
that 

(3.23) - — r — 27rlog —5 converges in <S'(R 2 ) as v — > 0. 

Moreover, away from the origin the convergence is as a symbol of order —2. Con- 
jugating with the Fourier transform yields an operator M(y) = (A R 2 + v 2 )^ 1 — 
l/27rlog( 1/^)1 whose kernel has a pseudodifferential singularity (that is, conormal) 
at the diagonal, but with growth as \z — z'\ — > oo. (Here, 1 denotes the operator 
whose kernel is identically equal to one, not the identity operator.) Let p be a 
function of compact support on R 2 , which is identically equal to one on a ball of 



large radius R. Then, since ( 3.23 ) converges as a symbol of order —2 away from 
the origin as v — > 0, pM{y)p converges as v — > as a pseudodifferential operator of 
order —2, and therefore, as a bounded map from H~ 1 (M. 2 ) to i/ 1 (R 2 ). 
Therefore, 

(3.24) R^ 2 )- 1 =T dirP (A R 2 +j, 2 )- 1 /9t /= ^log-l + T dir/ oM(i/)pJ, 

where 1 now denotes the operator on L 2 (H) with kernel equal to one. Because J 
is a continuous map from L 2 (H) to iJ~ x (R 2 ) and T^w is a continuous map from 
ff^R 2 ) to L 2 (H), the second term is a family of operators on L 2 (H) with uniformly 
bounded norm as v — > 0. Denoting the length of if by L, the operator 1 is L times 
a rank one projection on L 2 (H), so we see that the operator norm of i?(i/ 2 ) -1 is 
equal to 

n^ 2 r 1 n = ^iog- + o(i) 

as v — >• 0. Taking the logarithm of this, we see that 

log(A (A*)) _1 = log log/ 
Thus, logdeti?(/i) has the expansion 

(3.25) logdet-R(^) = - loglog/^ 1 / 2 -log^- + logdet 'i?(0) + o(l), /i 0, 

2n 



logCAoM)" 1 =loglog/i- 1 / 2 + log^+ (l), m^0. 



which is (3.5). This completes the proof of Theorem |2.6| . 

4. Compactness of Isophasal sets 
we defined the C°° topology on exterior domains. In this section 



In section 1.4 
we will prove 

Theorem 4.1. Each class of isophasal planar domains is sequentially compact in 
the C°° -topology. 

4.1. Compactification of the problem. First we explain how the problem is 
equivalent to a problem about a bounded, non-flat domain on S 2 . Consider an 
exterior domain Q whose closure does not contain the origin. Let g be a metric on 
the plane of the form 

/(M 2 ) 

where go is the flat metric, and f(t) is equal to one for t < R and f(t) = t 2 for 
t > R' . Then, g is the same as the standard metric on the ball B(R,0), whilst 
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its behaviour at infinity means that g extends to a smooth metric on S 2 regarded 
as the one-point compactification of R 2 . This compactifies ft to a domain ft' in 
S 2 . Let G be the composition of inversion in the unit disc of R 2 , followed by the 
identification above to S 2 . G is then a conformal map from R 2 to S 2 which maps 
the inversion fl 1 of £1 to ft'. Thus the metric g on ST pulls back to the metric 
e 20 o (w) d w( $yj on Q 1 5 w here w is a complex variable acting as a coordinate on R 2 in 
the usual way. Here = log \ G'\; note that </>o is not harmonic, because W is not 
flat. Let F be a conformal map from the unit disc D to ft 1 . Then the metric on 
Q' pulls back to e 2 ^ ' F (*" e 2 ^dzdz, where ^ = log is a harmonic function (see 



17 1, section 1). 

Now suppose that varies within an isophasal class. Since the topology on 
exterior domains is specified in terms of their inversions, the first thing we need to 
do is position Q well with respect to the inversion map. Recall that the first two 



heat invariants tell us that the perimeter and area of Q are fixed. By Lemma of 
the appendix, then, there is a uniform upper bound on the diameter, and a uniform 
lower bound on the inradius, over the isophasal class. Therefore, there is some 
r > such that every £1 in the isophasal class can be moved by a Euclidean motion 
so that the boundary lies in the annulus 

A r = {z e M 2 I r < \z\ < -}. 

r 

Consequently, the boundary of the inversion also lies inside A r . We choose R, in 
the definition of the metric g above, to be larger than 1/r so that all our domains 
fl can be placed isometrically within the flat part of the metric g. 

Next we compare our surgery formula for the exterior log determinant, 

(4.1) logdet'Ao -flog det A c + logdet R = 7 + log - 
with BFK's formula for H C S 2 : 

(4.2) logdet An' + log det A a + log det R = logdet A (S 2 g) +log^. 

Notice that the two R operators are the same, since the Laplacian is conformally 
equivariant and thus the harmonic extension of a function on H to the exterior 
is the same whether we use the flat metric or the metric g on the exterior of the 
obstacle. The quantities logdet A (52 and A are constant since we have fixed g 
once and for all, and as £1 varies over an isophasal class, the log determinant of An 
is fixed. Thus, subtracting the two equations we get 

log det An' = constant 

over any isophasal class. In addition, since the metric on S 2 is flat in a neighbour- 
hood of the complement of fi', the domains f2' all have the same heat invariants. 

Thus, our situation is that we have a class of metrics e 2 ^ ^^ e 2< ^ dzdz on the 
unit disc, which have fixed determinant and heat invariants. This is the same 
information as in flTj , but in this case the metric has an extra factor of e 2 ^ , where 
4>o is evaluated at the variable point F(z). In the next two subsections we adapt 
the argument of OPS to show that the set of <p's are compact in C°°, and in the 
final subsection we use this to prove sequential compactness. 

Remark. We can express the log determinant of An differently by considering 
also BFK's surgery formula for S* 2 with metric g with respect to the unit disc, D. 
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If we write D' for the complement of the unit disc, then this takes the form 

(4.3) logdet A D +logdet A( D >, g ) +logdeti?£> = logdet A.(s*,g) + 

Noting that logdet A^> and logdeti?u are universal constants, we get by adding 



(lyl) and ( |4.3| ) and subtracting ( |4.2| ) that 

logdet 'An = logdet An' — logdet Am', g ) + universal constants . 

Thus, up to universal constants, we can compute the log determinant of an exterior 
domain Q as a difference of log determinants on two bounded domains. One is for 
the domain Q' obtained by putting a metric on the plane that is Euclidean on some 
large ball and conformally compactifies the plane at infinity, and the other is the 
exterior of the unit disc with respect to the same metric. It seems likely that one 
could use this formula, together with a limiting process where g becomes Euclidean 
on larger and larger balls, to find an explicit Polyakov- Alvarez type formula for the 
exterior determinant. 

4.2. The Sobolev i estimate. Let us begin by recalling the way in which OPS 
proved C°° compactness for isospectral planar domains. Proving compactness is 
equivalent to obtaining uniform bounds on all Sobolev norms of the function <fi on 
the boundary of the disc D, which determines the isometry class of the domain O 
as the image of a conformal map F. Specifically, <f> extends to a harmonic function 
on the disc, and then F is the unique analytic function satisfying log = 4>, and 
-F(O) = 0, -F'(O) > 0. Once uniform bounds on the first few Sobolev norms of <f> 
have been obtained, it is easy to use Melrose's formulae for the heat invariants to 
obtain uniform bounds on the higher Sobolev norms inductively. 

As in OPS, we need to place an additional constraint on <p, namely that <f> is 
'balanced'. Since F is only determined up to a Mobius transformation of the circle, 
F, and therefore 0, are not uniquely determined. OPS made the definition that <$> 
is balanced if it satisfies 

e' p e l9 d6 = 0. 



s 1 

They proved that there is always a balanced conformal map from any domain, 
flat or not, to the disc |17| . This condition is important since then <fi satisfies an 
improved inequality, as discussed below. 

The crucial estimate of OPS is using the log determinant to obtain an if 1 / 2 
estimate on <fi. This goes as follows: using the constancy of the log determinant on 
the isospectral class, and the Polyakov- Alvarez formula for the log determinant in 
terms of <fi, we obtain 

(4.4) ~ f 09„0+ f <f> = C, 

z Js 1 Js 1 

Here and below, C, C±, etc, will denote constants which are uniform over the 
isospectral class. The first term is almost equal to the square of the Sobolev 1/2- 
norm. In fact, expanding <f> in a Fourier series, 



it is easy to show that 



(27T)- 1 / cpd n <j>d6 = Y J \n\\a n \ 2 , 
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SO 



d6+\a 2 \ 



[ <j)d n (j) de + 


[ <i>d6 


Is 1 


Js 1 



is equivalent to the square of the Sobolev i-norm. For future reference we note 

||^9„0|| i2 + 



that 
(4.5) 

is equivalent to the Sobolev (j + l)-norin. 

On the other hand, the Lebedev-Milin inequality for balanced <f> (Hii equation 
(5) of the introduction) gives 

1 



(4.6) 



logL 



< 



Combining the two we find a bound on J sl 4>d n (f> and then on | f sl <j>\, yielding a 
H 1 / 2 bound. (Without the balanced hypothesis (4.6) is only valid with coefficient 
\ in front of the final term, which does not yield any Sobolev bound.) 

In our situation, the metric e 2 ^ +2 ^° oF dzdz is not flat, so we get additional terms 
in our expression for the log det. Let us write <p t for <j> + 4>o ° F. As above, we may 
assume that <f> is balanced. Also, since the area and length are isophasal invariants, 
we have 



(4.7) 



e 0t = L, 



s 1 



= A' 



D 



where A' = Area(5 2 , g) — A and A is the common volume of the obstacles in our 
isophasal set. 

Lemma 4.2. There is a uniform bound on the Sobolev half-norm of both <fi and <pt, 
regarded as functions on S 1 , as fl ranges over an isophasal class. 

Proof. Constancy of the log determinant of fi' over the isophasal class implies that 
(111, equations (1.15), (1.16)) 



(4.8) / |V0 t | 2 

Expanding this out, we get 

|V0 + V(0 o oF)| 2 + 2 



t + 3 / d n <j> t = C 



oF + 3 / d n (<f> Q oF) = C. 



Since \a + b\ 2 > 3/4|a| 2 - 3|6| 2 , we have 
3 



D 



\V{<t> Q oF)\ 2 +2 + 2/ + 3 d n {<p oF)<C 



Integrating the first term by parts, using sup bounds on cf>o and \7(f>o, and using 
\F'\ = we get 

3 



3d 



D 



Since we have ( |4.7| ), and 



o 20 



,20 



2C 2 + 3C 3 



< C. 



e 20 t e -2<j, oF 



D 
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< c. 



the left hand side is bounded by e 2 "^ "°°A'. The term J gl is bounded similarly. 
Thus we have 

(4.9) ^ / <pd n cj> + 2 / 

4 Jd Js 1 

Adding twice ( |4.6| ) to this inequality gives us a bound on J <fid n (f>, and then (4.6) 
provides a bound on | J <j>\. This gives a bound on the Sobolev one-half norm of <f>. 

Next we bound the H 1 / 2 norm of <f> t . Since 4> t = 4> + 4>q o F, it is sufficient to 
bound the H 1 norm of tfio o F. The derivative of <f>o °F is bounded by \<f/ o F||.F'|. 
We have uniform sup bounds on <f>' . By Trudinger's inequality, (Pq[, (3.63)), 

(4.10) / e*<C||^|| ffl/2 , 

J s 1 

so 

lll^llla < hm<c\\<j>\\ 2 H1/2 <C, 

which gives us the required bound. □ 

4.3. Uniform L°° and Sobolev 1-bounds on <j). To obtain a uniform bound on 
1 4>\, we use the heat invariant a\, which shows that Jl7| 

k 2 (s)ds = / e-^{l + d n <p t ) 2 d9 = C. 

H Js 1 



Using a < 2(1 + a) + 2, we obtain from this 



e-^idnfafdO < 2C + 2 / , 
s 1 Js 1 



But we have a uniform bound on ||</>t||iji/2 so the Trudinger inequality gives a 
uniform bound on Je^^. Thus, 



e- 0t (5„0 t ) 2 ^<Ci- 

s 1 

2 /o _ f.2 ^ {„ I l\2 



Writing this in terms of </> and </>o, and using a /2 — 6 < (a + 6) , we obtain 



e 



({d^) 2 -2{d n (fr ) oF)) 2 )d6<C l . 



Since = e*, we get 



e-*(flUr <Ci + 2 e*((dn(h) ° F) z < C 2 , 
s 1 Js 1 

applying the Trudinger inequality again. The argument of p7[ , equation (1.13) 

on) can then be applied verbatim to conclude that sup</> and are uniformly 

bounded. 

4.4. Higher Sobolev bounds. Here we will prove that for any k, cf> is uniformly 
bounded in the Sobolev /c-norm. The proof is by induction on k. First we recall the 
proof in the OPS case, where we just have the harmonic function (f>. Thus suppose 
that H^Hffj is uniformly bounded, and, therefore, H^Hc^- 1 is also, by the Sobolev 
inequality. We start at j = 1 since the hypotheses have been proved for this value 
above. 

We use Melrose's result that 

' (dik) 2 (s)ds 

H 
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is uniformly bounded. This implies that 

/ {dlkfdO 
Js 1 

is uniformly bounded, since 

(i)' - (•-£)' (!)'+|><«-'.**-.*- i <, Vd9 

where pi are polynomials, and all the occurrences of (f> above are -L°°-bounded by 
the inductive assumption. 
Consider 

in view of (4.5), it is sufficient to bound the L 2 norm of this quantity. Taking 



derivatives, we find that this is equal to 

(4.11) (d^)k + Y^ Cl (dy)(df l k). 

1=0 

The term 9™e^ is equal to 

(8?(f>)e+ + q(e+,de(i>,...,82 > - 1 4>), 
with q a polynomial, and is therefore uniformly L 2 -boundcd for m = j, and uni- 



formly L°°-bounded for m < j. Thus, the L 2 norm of (4.11) is bounded by 

C(||^e*||2||*|| a + X;il^ll~ll^" , *lla)' 
which is a uniform bound. 

In the exterior domain case, uniform Sobolev bounds on <j) t follow in the same 
way. Now we obtain uniform Sobolev bounds on <j>. Since 

(4.12) %d n <t>t = 4 d n<f> + didnifo o F), 

we need to obtain a uniform L 2 -bound on the second term. Recall that this depends 
on cf> through the function F, since <f> — log \F'\. Notice that 

d n F = e l6 , 

where tp is the harmonic conjugate of ip. Thus 

4d n (^ o F) = di((d n ct> ) o F{e w e^)) 

which involves at most j derivatives of <p and ip. Since tp is the harmonic conjugate 
of cf>, normal, resp. tangential derivatives of ip are equal to tangential resp. normal 
derivatives of 4> up to factors of i, and can therefore be estimated by derivatives of 
ip. The only way that j derivatives of <fi or tp can occur is in the term 

({dnMoF)(dl{e w e^)) 

which can be estimated by 1 1 c?ti <^»o 1 1 oo | i f?g lis- The other terms can be estimated 
by II^oIIcj 'll^e" ^llci- 1 ■ This shows that the second term of (4.12), and therefore 



also the first, is uniformly bounded in L 2 . This completes the inductive step of the 
proof. Thus, (j> is uniformly bounded in C k , for all k. 
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4.5. Proof of Theorem |4.l| . Consider any sequence of fi^ of isophasal exterior 
domains. We will show that there is a subsequence converging to some domain f2 
in the same isophasal class. 

First we show that there is a convergent subsequence. As discussed above, we 
may assume that each Qk is placed so that its boundary lies inside some annulus A r . 
Let fi£ be the inversion of Q^. Then we have shown above that any set of balanced 
4>k corresponding to the f2£ is precompact in the C°° topology. Consequently the 
set of normalized conformal maps F k (given by <fik = log \Fj,\, F(0) = 0, F'(Q) > 0) 



is precompact in the C°° topology |17 . Passing to a subsequence, we may assume 
that the Fk converge in C°° to a conformal map F. 

The domain f2£ is then given by f2£ = Ek(Fk(D)) for some Euclidean motion 
Ek- Note that all £l£ lie inside the ball of radius R. Thus the collection {Ek} of 
isometries lie inside some compact set; in fact, since Ek(0) is contained in 5^(0), 
we have 

E h c{T i oR \6e[O,2ir],\b\<R} Vfc 

where Tf, is translation by complex number b and Rg is rotation about the origin 
by angle 8. Therefore, some subsequence of Ek, say Ej h , converge to a limiting 
Euclidean motion E. Along this subsequence, Fk = Ej k o Fj k converges in C°°, so 
the sequence of domains ilj k converges to SV = E(F(D)) . Thus, by definition of 
the topology on exterior domains, the sequence ilj k converges to f2, the inversion 
of H, 1 . Thus the class of isophasal domains is precompact. 

To complete the proof we have to show that the isophasal class is closed in the 
C°° topology; that is, if a sequence of obstacles ilk with the same scattering phase 
converges in C°°, then the limiting obstacle fl has the same scattering phase. 

The scattering phase is given by 

s(A) = -i log det S(X) = -Oogdet (ld + J^A(X)), 

V Z7T 

where A(X) is the operator on L 2 (5' 1 ) with C°° kernel 

(4.13) A{\)(6, U ) = I ue^^ - 

J\z\=r ov dv 

for sufficiently large R (see the appendix of Q). Here uq is the distorted plane 
wave with incoming direction 8. Let us fix A. The kernel ^4(A) belongs to the trace 
class of operators and s(A) is a continuous function of ^4(A) with respect to the 
trace norm. Thus, we need to show that ^4(A) is a continuous function, as a trace 
class operator, of the domain. It suffices to show that the distorted plane waves ug, 
together with a certain number of derivatives in 8, are continuous functions of the 
domains, say in Lf oc . We will just treat the case of ug itself, since the ^-derivatives 
are treated similarly. 

Choose a C°° mapping Ik : Q — * Clk which is the identity for \z\ > 2R. This can 
be done so that the maps Ik converge in C°° to the identity map. Pulling back the 
operator A — A 2 to f2 gives us differential operators Lk — A 2 on fi, with Lk = A 
for | z | > 2R, and with coefficients converging to those of the Laplacian in C°° as 
k — > oo. Thus, moving to S7, we must find solutions to 

{L k - \ 2 )u k =0, lit = on u k - e~ lXz ' e outgoing 
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and show that u k —* u, where u solves the limiting problem on f2. By standard 
methods this reduces to solving 

(Lk - X 2 )u k = g k , u k = on dQ, u k outgoing 

with g k — > g in C°°, and proving that u k — > u, where u solves the limiting problem. 
These equations may be solved using the method of the appendix of Q (which in 
turn is based on the proof of Theorem 5.2 of p|). Tracing through the proof, it is 
not hard to show that indeed u k — > u in Ll ' . 

Therefore, the scattering phase s(A) for Q, k converges pointwise to the scattering 
phase for f2. This shows that fl is isophasal with the £l k and completes the proof 
of the theorem. 



5. Appendix 

5.1. Uniform bounds on isophasal classes of domains. In the proof of The- 



orem 4.1 we need the following lemma; recall that any class of isophasal domains 



has fixed area and perimeter. 

Lemma 5.1. For any obstacle O there is a lower bound on the inradius, and an 
upper bound on the circumradius, depending only upon the area and perimeter of 
O. 

Proof. The bound on the circumradius R < L, where L is the perimeter, is trivial. 
Consider the inradius. Let O be a domain with area A and perimeter L and r a 
number such that r/2 is larger than the inradius but smaller than the circumradius. 
It is possible to choose of covering of O by balls of radius 5r whose centres lie in O, 
such that the balls of radius r with the same centres are disjoint. To construct such 
a covering, start with a finite covering by balls of radius 3r whose centres lie in O. 
Choose any two of the balls. If the balls of radius r about their centres intersect, 
then the ball of radius 5r about any one of the centres contains both balls of radius 
3r, and therefore one of the balls can be discarded. By discarding balls successively 
in this way, we end up with a covering with the required property. Let N be the 
number of balls in the covering. 

Considering the area of each ball, we have an inequality 

(5.1) 257rr 2 iV > A. 

Since r/2 is assumed larger than the inradius, but smaller than the circumradius, 
the circles of radius r/2 with the same centres as the balls in our covering must all 
intersect the boundary. Since the circles are all distance at least r apart, this gives 
an inequality 

(5.2) L > Nr. 
Combining the two inequalities, we find that 

25?rr > — . 
L 

Taking the infimum over r, we find that 

2A 



inradius(O) > 



25ttL 

□ 
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5.2. Computation of the scattering phase for the unit disc. To compute 
the scattering phase, we construct the distorted plane waves for small A. These are 
given by 

u u , x (z) = e~ iX *-" + < A (z), zeR 2 , ujeS 1 
where u' x satisfies 

(A + A 2 )< A = 0, 
(5.3) u'^ x {z) \ S 1 = -e~ iXz -« r S 1 , 

< <x {z) = M^V^^A, z, -co) + 0{\z\^l 2 ), \z\ -> oo. 

Then the scattering phase is equal to — i logdet 5(A), where the scattering matrix 
is the unitary operator 

5(A) = Id+ J^A(X), 



and A(X) has kernel (obtained by applying stationary phase to (4.13)) 

A(X)(9 1 ,9 2 )^e l ^ 4 a(X,9 1 ,9 2 ). 

To find u'^ A , we decompose e ~ tXz ' u \ S 1 as a Fourier series; u'^ A then has a 
corresponding decomposition into Bessel functions. 

Let 9 = B\ — 02, where 9i is now regarded as a circular variable in [0, 2tt). By 
circular symmetry, A(X) depends only on 9. In terms of 9, on the unit circle we 
have 

— iXz-uj iXcosO 

e = e 



(ixy(e w + e - w y 

Hence, on the unit circle, 



^ 2'/! 

i=o ■' 



^—iXz-uj 



= J2 a n {X)e me , where 



71— — OO 



(5-4) a Q (X) = 1 + A 2 a (A), a (A) bounded, |A| < 1, 

K(A)| < ™ |A| < 1. 



n ! 



Then 



Va( z ) = 2^ fl n(A) ^ e ,0 = z-(-w), 

-n I ^ I V / 



where if„ is the Hankel function of the first kind of order n. Using large \z\ 
asymptotics of the Hankel function ((!]], chapter 9), this gives 

, — in/ 4 — in.Tr/ 2 in9 



ttA J?|„|(A) 



Since |i? n (A)| < A" for A < 1, and we have the estimate (5.4) for a„(A), this series 
converges. As A — * 0, 

# (A) = -logA + O(l), 

7T 
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so we obtain for a(A, 9) the asymptotics 



and for A(X), 



a(X, 9) = J^^ 4 ilg A + 0(X~ x ^ 2 (ilg A) 2 



A(A)(0) = -ilgA + O((ilgA) 2 ), 



where the error term is a smooth function of 9. Thus, 

s(A) = -i logdet 5(A) = -i tr log ( Id +A{\)) = it ilg A + 0((ilg A) 2 ). 



5.3. Proof of Lemma |3.l| . The operator (A ffi + /i) 1 — (A R 2 + /i) 1 is given in 
terms of the heat kernel by 

/•oo 



(5.5) / e-^( e - tA ® -e~ tA ^)dt. 

Jo 

Jensen and Kato proved the following estimate of the trace norm of the difference 
of heat kernels: 

||( e -*Ae _ e— * A »»)l|x =0{t- 1 ' 2 ), t->0. 
On the other hand, if we denote the two semigroups by H(t) and Ho(t), then 

(5.6) \\H(2t) - flo(2t)||i = \\H(t)(H(t) - tf (*)) + (H(t) - if (*)) #o(*)l|i 

< \\H(t)\\ op \\H(t) - HoWlli + II^W - ffo(*)||i||#o(*)||op = 2||H(t) - HoWlk- 

Iterating this ineq ualit y shows that the trace of H(t) — H (t) is O(t) as t — > oo. 
Thus, the integral (5.5) is convergent in trace norm, and so the result is trace class. 

Next we prove the second part of the lemma. Using the functional calculus, we 
have 

( e -T(A ffl +p) _ e -r(A B2+At )) dr = + ^_l e _t(A e +,.) _ ( Ar2 + fl yl e -t(A R2+f .)_ 

Using the estimates on trace norms above, we see that the right hand side is trace 
class, the trace is differentiable as a function of t, and the derivative is minus the 
integrand on the left hand side of the equation. Hence we can calculate, for 5Rs > 1, 



(5.7) 



dfi ' ' dfi r(s) jo ^ 

1 Z" 00 

~ / t s tr ( e -*( A e+f) _ 6 -*(A e2 +m)) dt 

r(s) jo 



i f 00 ^ 



r(s) ^ ■ |tr((A e + M )- 1 e ^ A ^)-(A R2 + A1 )- 1 e- t ( A « 2+ ^)^ 

= ~ ffe / °° <S tr ( (Ae + A*)-^-*^ 8 ^ - (A R2 + M )- 1 e-*(^+")) |. 

The boundary term in the integration- by-parts is zero since, for 5Rs > 1, the inte- 
grand tends to zero at both zero and infinity. Since (A® + /i) _1 — (Ar2 + /i) _1 is 
trace class, it follows that the integral has a simple pole at s = 0. But the term at 
the front, s/T(s), has a double zero at s = 0, so when we take the minus derivative 
at s = to find the derivative of the log determinant, we obtain precisely the pole of 
the integral at s = 0. This pole is tr((A® + /i) _1 — (A R 2 + /i) -1 ), proving (3.6). □ 
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5.4. Hadamard-regularized integrals and Pushforward. Let h(x) be smooth 
and compactly supported on [0, oo). Then J hdx/x is not convergent. We define 
the Hadamard-regularized integral of h by the limit 

f°° dx / f°° dx 1 

HR-/ h(x) — = lim ( / h(x) — - h(0) log - 
Jo x £^oV/ £ x e 

It is easy to check that the limit exists. It may also be described as the constant 
term in the asymptotic expansion of J hdx/x as e — > 0. 

Hadamard-regularized integrals turn up naturally in the Pushforward theorem 
for poly homogeneous functions proved by Melrose |L3|, ||. In fact, the authors 
first derived the expansions in section 3.3 using a special case of this theorem, so 
we will include a brief discussion. 

The pushforward is invariantly defined on densities rather than functions, so we 
consider densities defined on R^_. It is most natural to consider b-densities, that is, 
densities of the form 

. . dxi dx2 

g(xi,x 2 ) . 

Xi x 2 

One reason for this is that such densities have an invariantly defined restriction 
to the boundary faces Xi — 0, obtained by cancelling the factor dxi/xi (which 
is invariant under changes of boundary defining function at Xi — 0). Thus, the 
restriction of g(xi,x 2 ) dxi/xi dx 2 /x 2 to x\ — is g(0, x 2 )dx 2 /x 2 . We will consider 
g which are smooth and have compact support. Then we have 

Proposition 5.2. Consider the map f : — > R+ given by 

{xi,X 2 ) l-> X = X\Xi- 

Let u = v(xi,x 2 ) dxi/xi dx 2 /x 2 be a smooth b-density with compact support. Then 
the pushforward of u has an asymptotic expansion 

oo 1 dec 

(5.8) /*u~^( Pi +qr,.log-)sE J — , 

j=o X X 

where q = v(0, 0) and 

poo j poo 1 

(5.9) po=(HR-/ v(x u 0)— +HR-/ v(0,x 2 )—). 

Jo x\ J x 2 

This theorem is a special case of the general result in Jl3| , and is proved explicitly 
in [§, section 2. The expansion ( |3.22| ) follows immediately from this theorem by 



regarding the integral as a pushforward under the map (a, /3) i— > r = a(3 (we also 
have to multiply by the formal factor dr/r to make the integrand into a density). 
In fact, the expansion ( 3.18] ) can also be deduced from the pushforward theorem 



by using the operations of logarithmic and total blowup discussed in M . 
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